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Introduction 

At  the  Air  Force  Weapons  Laboratory,  Interest  has 
continued  for  some  time  in  energetic  plasma  injectors. 
A  possible  scheme  for  such  a  device  is  the  plasma 
deflagration  gun.  When  the  question  arose  whether  it 
would  be  possible  to  scale  a  deflagration  gun  to  the 
multi-mega joule  energy  level,  it  became  clear  that  a 
scaling  law  which  described  the  gun  as  a  circuit 
element  and  allowed  one  to  confidently  scale  gun 
paramenters  would  be  required.  We  sought  to  develop  a 
scaling  law  which  self-consistent ly  described  the 
current,  magnetic  field,  and  velocity  profiles  in  the 
gun.  We  based  this  scaling  law  on  plasma  parameters 
exclusively,  abandoning  the  fluid  approach  used  by 
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energy)  B*  =  0  at  some  station  of  the  gun.  Further, 
we  could  simply  extend  the  conductor  such  that  Er  =  E 
where  B  *  =  0  and  this  mechanism  would  imply  an 
infinite  velocity  for  the  plasma.  Since  this  is 
unphysical,  v  must  plateau  at  some  value  v  =  vf;o„,  .  A 
possible  explanation  of  this  cutoff  in  velocity  is  the 
particles’  Larmor  radius  exceeding  the  interelectrode 
spacing  (rL  >D/2),  but  for  purpose  of  this  discussion, 
it  is  sufficient  to  say  that  some  mechanism  does 
exist. 

We  set  z  =  0  at  the  point  where  acceleration 
ceases  and  where  v  =  vf;nat  .  We  identify  the  magnetic 
field  at  that  point  as  B,  r  B„,  so  that 


Scaling  Law 

Our  one -dimensional  quasi-static  plasma 
deflagration  gun  is  shown  in  Figure  1.  The  boundries 
at  y  =  0  and  at  y  =  D  are  perfect  conductors.  The 
boundry  at  y  =  0  is  maintained  at  a  constant  voltage  V 
with  respect  to  the  opposite  boundry  so  that  a 
constant  electric  field  Ey  exists  throughout  the 
volume  of  the  gun.  The  flow  of  current  through  the 
yet  unspecified  plasma  within  the  volume  of  the  gun 
yields  a  current  density  Jy(z)  and  a  magnetic  field 
B«(z)  within  the  volume.  The  current  path  enters  at 
y  =  0,  z  =  zb,  (zb>0),  crosses  between  the  electrodes 
throughout  the  volume  of  the  gun,  and  exits  at  y  =  D, 
z  =  zb  so  that  the  magnetic  field  is  zero  at  the 
muzzle  of  the  gun,  z  =  zm,  (zm<0)  and  the  plasma 
moves  in  the  -z  direction. 

We  begin  by  assuming  a  low  density  flow, 
dominated  by  single  particle^dynamios  in  which  plasma 
motion  is  driven  by  the  E  x  B  drift  of  the  particles 
as  they  orbit  the  magnetic  field  lines.  The  bulk 
plasma  velocity  is  given  by 


Now,  unless  the  current  distribution  extends  to 
z  =  -  oo  (which  would  require  infinite  magnetic  field 


From  z  =  0  to  the  muzzle  of  the  gun  at  z  =  z„,  the 
only  current  flowing  is  the  conduction  current 
J  =  o*bE  resulting  from  the  fact  that  the  plasma  has 
some  finite  temperature  and  hence  some  finite 
conductivity  er*v,  .  We  assume  that  the  plasma 
temperature  and  conductivity  are  constant  for  all  z 
(including  z  >  0)  so  that  Jc  is  constant.  For 
completeness  we^assume  that  the  power  input  into  the 
volume  by  Jc  x  E  is  balanced  through  radiative  and 
conductive  losses.  This  is  necessary  if  we  are  to 
find  a  steady-state  solution  to  the  problem,  and  is 
consistent  with  the  constant  temperature  assumption. 

If  particles  enter  the  gun  at  positive  z  and  are 
accelerated  until  they  reach  the  cutoff  point,  there 
must  be  an  additional  input  of  electrical  power  into 
the  volume  for  z  <  0.  We  allow  this  additional  power, 
input  through  the  presence  of  a  current  density  Jj 
such  that  Ja  x  E  contributes  only  to  the  drift  of  the 
particles,  and  it  is  this  current  density  for  which  we 
wish  to  find  a  self-oonsistent  solution.  We  begin  by 
describing  the  dependence  of  B(z)  on 
J(z)  =  Jc  +  Jj(z).  We  will  solve  the  problem 
quasi -statically  and  one-dimensionally ;  that  is,  we 
will  allow  variations  on  the  z-direction  only, 

S  1  37  =  St  1  0  (3) 
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Figure  1.  1-dimensional  quasi-static  plasma 
deflagration  gun. 


It  is  easy 
B(z)  = 


so  that  we 

dB(z) 

dz 


to  see  from  Ampere's  law  that  for  z  >  0, 

B.  +  m.JcZ  +  Jd(z')dz'  (4) 

Jo 

may  write 

=  +  m0J«j(z)  (5) 


Now.^if  we  say  that  the  bulk  plasma  velocity  is  simply 
the  E  x  B  drift  velocity  at  any  station  in  the  gun,  we 
have 


v2(z) 


E* 

FTzT 


and 


dv*(z) 

dz 


..  Ea  dB(z) 

B*77T  dz 


(6) 


(7) 
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We  choose  to  evaluate  v!  because  it  is  related  to  the 
directed  kinetic  energy  of  the  ions,  by  which  the  bulk 
of  the  momentum  and  energy  are  carried.  Per  particle, 


*<*>  ■  (f 

L)  v  2 (z ) 

(8) 

dK(z)  /m 

;\  dv2(z) 

(9) 

dz  '  \2 

1  dz 

Now,  if  the  particles  are  streaming  through  the  gun  so 
that  the  particle  flux  is  4>,  the  differential 
electrical  power  per  unit  length  input  into  bulk 
plasma  motion  (recalling  that  particle  energy  is 
increasing  with  decreasing  z)  is 


dP(z ) 
dz 


=  -E  Jj (z) A 


(10) 


But  the  differential  power  input  must  be  related  to 
particle  kine 

„dK(z) 


4-A¬ 


dz 


energy 

by 

«(?■ 

\  dv2(z) 

1  dz 

(id 

Mr! 

dB(z) 

dz 
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| 
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|  Jj  (z '  )dz| 

'  j^MoJc  +  MqOj  ( z  )j 
Equating  the  two  expressions  yields 


03) 


-E  Jd(z)A  =  -4>Am 


;E2|e 
**oJc  +  VoJ4  (z 


B0  +  *■  U 


Jd(z')dz'j 


(14) 


Solving  for  J  (z),  we  arrive  at  the  integral  equation 


*Jj  (z)  =  Jc 


m;4>E/u0 


rz  13 

B0+  n0Jc  z  +  Mo\  jj(z')dz  -  m;4>EM0 

L  Jo 


(15) 


This  equation  may  be  evaluated  numerically  to  arrive 
at  a  current  profile,  but  we  may  readily  solve  for  J 
at  z  =  0  where  acceleration  ceases. 


J,(0) 


.  m;4>EMo 
*  Ba  -  m(4>EMo 


(16) 


Note  at  Jj(0)  becomes  quite  small  for  small  m;<tE^0, 
and  becomes  negative  (corresponding  to  energy 
extracted  from  the  plasma,  an  unphysical  solution)  for 
m;4>E^0  >  B|  .  Therefore,  we  must  operate  the  gun  such 
that 


mi4>E„0  <  B3 


(17) 


Recalling  that  v^ma|  =  E  /B0 ,  we  rearrange  terms  to  get 


m;4-v+;Ml  < 


08) 


In  other  words,  the  pressure  (or  momentum  flux)  of  the 
plasma  stream  cannot  exceed  twice  the  magnetic 
pressure  at  cutoff.  It  is  noteworthy  that  this  is 
similar  (within  a  factor  of  two)  to  the  result  whioh 


Cheng  derrived  by  eliminating  terms  of  the  momentum 
equation  by  estimating  orders  of  magnitude.2  We  may 
define  an  effective  beta 


= 


m;4>v-,„'»t 

(B0V2^0) 


09) 


For  a  deflagration  gun  operating  with  maximum  output 
of  momentum  flux  (maximum  effective  beta),  Jd(0)  is 
infinite;  that  is,  power  input  into  the  bulk  plasma 
motion  takes  place  in  an  infinitely  thin  current 
sheath.  For  finite  Jj(0),  the  current  sheet  has 
finite  thickness,  but  an  examination  of  equation  (18) 
shows  that  the  larger  Jj(0),  the  quicker  J<j(z)  falls 
off  with  increasing  z.  In  short,  our  simple  model 
returns  a  result  that  has  been  verified  in  experiment: 
the  current  distribution  in  a  deflagration  gun  is 
composed  of  a  diffuse  current  throughout  the  barrel 
and  a  thin  current  sheath  of  high  current  density  at 
the  breech.  (See  Figure  2.) 


Comparison  with  Experiment 

Let  us  consider  the  B-field  profile  that  would  be 
observed  in  a  deflagration  gun  operating  at  maximum 
efficiency  (0e{f  =  2),  accelerating  particles  from 
rest  to  v  =  v*;„ai  in  a  infinitely  thin  jump  region. 
Granted,  this  limit  is  unattainable,  but  it  is  the 
simplest  to  analyze.  We  integrate  over  the 
delta-function  current  profile  Ja  to  obtain  the 
current  in  the  jump  region  Id.  Power  balance  requires 
that 

E  DIj  =  (f1)  v^.,  +  A  (20) 

Now,  we  may  use  our  two  equivalent  expressions  for 
Vf  ,  ,  equation  (2)  and  a  derivation  from  equation 
(19), 


along  with  Ampere's  law, 

AB  =  (22) 

A 


to  solve  for  the  step  increase  AB  across  the  current 
sheath.  Combining  equations  (2) , (20) , (21 ), and  (22) 
yields  the  result 

AB  =  B  (23) 
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Experimenter 
or  Group 

Charge 

Voltage 

Maximum 

Current 

Inner 

Radius 

Outer 

Rad ius 

Aspect 

Ratio 

Gas 

Species 

Output 

Velocity 

Output 

Density 

D.  Y.  Cheng 

20  kV 

0.65  MA 

0.63 

era 

4.5  cm 

7.6 

He 

3-12  omAus 

10  s  -  10" 

cm 

Case  Western  6 

9  kV 

- 

0.50 

era 

4.0  cm 

9.1 

He 

10-12  om/pis 

10's  -  10“ 

cm 

AFWL  5 

18  kV 

0.42  MA 

0.65 

cm 

2.5  cm 

12.3 

H,He,Ar 

11  cid/ms 

3-5  x  10,s 

cm 

Typical 

20  kV 

0.50  MA 

0.60 

cm 

4.0  cm 

10.0 

He 

1 1  cm/MS 

5  x  10,s 

cm 

Table  1.  Comparison  of  parameters  of  several  plasma  deflagration  guns. 
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Figure  5.  MAGPIE  1- dimensional  problem. 


Conclusion 


We  have  developed  a  simple  model  of  the  plasma 
deflagration  gun  which  includes  a  self-consistent 
description  of  magnetic  field  and  current  density 
within  the  volume  of  the  gun.  The  results  of  this 
analysis  are  in  agreement  with  data  available  from 
deflagration  gun  experiments.  The  analysis  points  to 
high-aspect  ratio  designs  for  deflagration  gun 
experiments  in  high-energy,  high-current  regimes. 
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